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We investigate phase coherent ballistic transport through antidot lattices in the generic case where 
the classical phase space has both regular and chaotic components. It is shown that the conductivity 
fluctuations have a non-Gaussian distribution, and that their moments have a power-law dependence 
on a semiclassical parameter, with fractional exponents. These exponents are obtained from bifur- 
cating periodic orbits in the semiclassical approximation. They are universal in situations where 
sufficiently long orbits contribute. 



I. INTRODUCTION 

Experiments on mesoscopic semiconductor devices 
have exhibited a variety of features that can be attributed 
to quantum chaosi^. Modern fabrication techniques 
make it possible to produce extremely pure semiconduc- 
tor microstructures in which the motion of the electrons 
is confined to two-dimensional domains whose relevant 
dimensions are much smaller than the phase coherence 
length and the transport mean- free-path. In this ballistic 
regime semiclassical methods have been very successful 
in connecting quantum interference effects to the under- 
lying classical dynamics 4 *^. 

There has been a considerable emphasis on systems 
in which the confinement potential leads to chaotic mo- 
tion. As in disordered systems, certain transport prop- 
erties are found to be universal, like conductance fluc- 
tuations and weak localization properties of transport 
through cavities§i2i2ii2iii. In contrast, if the dynamics 
is integrable, these features are in general not universal 
but depend on the specific system. Semiclassical meth- 
ods have been applied to explain universal properties of 
chaotic transport as well as the non-generic behaviour of 
integrable cavities&LSiiSiii 

Other experiments have been concerned with re- 
vealing signatures of classical periodic orbits in quan- 
tum phenomena. Examples are orbital magnetism in 
ballistic microstructures^ and transport through an- 
tidot superlattice3i£ii§iiLi2ii£. The latter consists of 
a two-dimensional electron gas at the interface of a 
GaAS/Al^Gai-^As heterostructure into which a peri- 
odic array of holes is drilled. The effective potential is a 
periodic structure of high potential peaks, and if the po- 
tential is steep it may be considered as an experimental 
realization of the Sinai billiard. Experiments on ballis- 
tic transport at temperature T ~ 1.5K show a series 
of pronounced peaks in the longitudinal conductivity vs 
magnetic field that can be explained by classical electron 
transport within the Drude formalismi£*i£. Experiments 
at a lower temperature T as OAK reveal additional quan- 
tum oscillations superimposed on the classical peaksi 7 -. 
These quantum oscillations can be attributed to unstable 
periodic orbits of the electrons in the confinement poten- 
tial by a semiclassical theory for the conductivitjsiSiiS. 



Generic systems have a phase space which is neither 
completely chaotic nor integrable but contains a mixture 
of regular islands and chaotic regions. This is relevant 
for experiments because confinement potentials are not 
hard-wall potentials, and soft-wall potentials typically 
lead to mixed dynamics. A natural question is whether 
ballistic mesoscopic systems have characteristic proper- 
ties that differ from those of chaotic and integrable sys- 
tems. One signature of mixed dynamics has been found 
in the ballistic transport through cavities. It was pre- 
dicted and observed experimentally that the variance of 
conductance fluctuations, as a function of the magnetic 
field, has a power law dependence with a non-integer ex- 
ponent that is related to the trapping of chaotic trajecto- 
ries near regular islands 20 - 21 . The purpose of this article 
is to show that there is a different mechanism in ballistic 
transport through antidot lattices which leads to univer- 
sal behaviour in mixed systems. 

Amongst the main characteristics of the dynamics in 
mixed systems are bifurcations of periodic orbits, events 
in which different periodic orbits coalesce when param- 
eters of the system are varied. Bifurcations are impor- 
tant for semiclassical approximations because bifurcating 
orbits carry a semiclassical weight that is higher than 
that of unstable periodic orbits and sometimes even that 
of tori of regular orbits. The dominating influence of 
bifurcations on transport through antidot lattices has 
been demonstrated in 2 ^ 3 .. Bifurcations occur in differ- 
ent forms, and depending on the physical quantity that 
one considers, they are of different importance. In the 
following we consider moments of conductivity fluctua- 
tions. We show that they are dominated by a competi- 
tion between different types of bifurcations, leading to a 
power-law dependence on a semiclassical parameter with 
fractional exponents. These exponents are universal if 
the competition is amongst all generic bifurcations. Sim- 
ilar results have been obtained for moments of spectral 
counting functions and wavefunctions in closed mixed 
systema 24 ! 25 ) 26 ! 27 . 

In the following we briefly review the semiclassical 
theory of transport through antidot lattices and discuss 
modifications in the presence of bifurcations. We give 
an overview of different types of bifurcations and discuss 
their influence on moments of conductivity fluctuations. 
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II. SEMICLASSICAL THEORY FOR THE 
CONDUCTIVITY 

The starting point for the semiclassical theory is the 
Kubo linear response theory. The conductivity is given 
in terms of matrix elements of the current operator. The 
semiclassical approximation of the matrix elements^ and 
the application of the stationary phase method for eval- 
uating integrals then yields a semiclassical expression 
for the conductivity as a sum over the classical (Drudc) 
component^ and an oscillatory component in terms of 
periodic orbit a 18 i 19 . 

The Kubo formula for the longitudinal conductivity 



~v 



Tr [v x 5 r (E F - H)v x 5 T {E F - H)} (1) 



where H = (l/2m)(p — eA(f)) 2 + U(r) is the Hamilto- 
nian of a two-dimensional electron gas in a perpendicular 
magnetic field, A = (— By/2, Bx/2) is the vector poten- 
tial taken in the symmetric gauge, U(r) is the confine- 
ment potential, and v x is the x component of the velocity 
operator. E F denotes the Fermi energy, m is the effec- 
tive mass of the electron, and V is the area of the system. 
Weak disorder is taken into account at the level of Born 
approximation by giving the <5-functions a finite width 
r = h/2r e i where r e j is the elastic scattering time. 

To evaluate the Kubo formula it is conve- 

nient to write the delta function in terms of re- 
tarded and advanced Green functions: 5r(E F — H) = 
-(l/2m)[G+(E F ) - G-(E F )\ with G ± (E) = [E - H ± 
£T] . In the semiclassical limit as h — > 0, the longitu- 
dinal conductivity is reduced to terms that involve only 
products of the retarded and advanced Green functions, 
because terms that involve products of Green functions 
of the same type vanishi^ 1 *^ 



e 2 Trh ( 1 



d 2 rd 2 r'v x (r)v x (r')x (2) 



V \2ir 
G- T ,{E F )G+ !T (E F ). 

The Green function G^ r , r s (E) is semiclassically approx- 
imated by a sum over classical trajectories 7 that start 
at the initial position r and end up at the final position 
r' with energy E given bj*2fl*2i: 



Gi r — 



exp 



2irh? 



2r, 



1 h ^ r ~ l 2^ 



(3) 



S-y(r' , r) = J pdr is the classical action of the trajectory 
7 with the canonical momentum p, rj~ t is the number of 
conjugate points along 7, T 7 is the traversal time of the 
trajectory and the first exponential is the weak disorder 
factor. The amplitude factor D 7 involves second deriva- 
tives of S 1 . It is convenient to express it in terms of local 



coordinates r = (z, y), with coordinate z along the orbit 
and y perpendicular to it. Then 



vv' dydy' 



(4) 



where v and v' are the velocities at z and z' . Inserting 
(J3J) into J5J one finds that the longitudinal conductivity 
is given by a double sum over oscillatory terms. Due to 
the fact that G + (E) and G~(E) are complex conjugates, 
G~ r ,(E) — [G*, T (E)]*, the phase of the oscillatory terms 
contains the action difference of two trajectories (71 and 
72). The conductivity is then split into two parts, a non- 
fluctuating one from the diagonal terms (71 = 72 = 7) 
and a fluctuating one from the non-diagonal terms (with 
7i 7^ 72), that is: 



^xx ~ ®xx ~r~ ^xx' 



(•5) 



The mean conductivity a xx 

e 2 f \ 

°** = h 2 V / d2rd2r ' v*(r)«x(r')-D 7 exp(-T 7 /r e O, 

7(r,r') 

can be transformed into a familiar form as a phase-space 
average (denoted by (-) r ,p) where it is recognized as the 
Drude conduct ivitjsi&iS 

/>oo 

0^ = e 2 p(E F ) / dt (v x (0)v x (t)} r , p exp (-t/ T e i) (7) 
Jo 

where p(E F ) = m/2-Kh 2 is the density of states per unit 
area at the Fermi energy in two dimensional systems. 

The fluctuating part a! xx is a quantum correction to 
the Drude conductivity and is given by the semiclassical 
formula 



71,72 



d 2 r d 2 r'{v l ) x {v' 2 ) x ^D~D 1 



T \ ( % Z7T 

' M> ! - — 277; exp 1 ^ s ^^ r > r '> EF "> ~ y^i,72 / 
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where T~, 



(T 7l + T 7l ); S. 



5 7l and 



7l ,72 

7 77i,72 — Vfi ~ V~t2- A correct evaluation of these in- 
tegrals requires a detailed knowledge of classical phase 
space structures as will become clear as we proceed. 

The main contribution to the integrals comes from the 
stationary points where 



V5 71)73 (r',r;E F ) = p 2 - Pi = 
V'5 71)73 (r' ! r;E F )=p' 1 -p' 1 =0 



(9) 



Here Pi and p' ' i denote the initial and final momenta 
of trajectory 7$. The two trajectories 71 and 72 must 
hence have the same initial and final momenta in addition 
to having the same initial and final positions. One way 
of satisfying these conditions with 71 and 72 not being 
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identical, is that both trajectories are part of a primitive 
periodic orbit 7. In fact, for each periodic orbit 7 there 
is an infinite set of pairs that satisfy the stationary-phase 
conditions. They differ only by the number of times they 
wind around 7 and their action difference is a multiple of 



etition of a bifurcating orbit 7 is 



the action of the primitive periodic orbit SU 



-S 7l — 



lS>y 



The evaluation of JHJ is closely related to performing 
the trace in the derivation of Gutzwiller's Trace Formula. 
The integration is done in the local coordinates r = (z,u) 



with d 



dzdz' dydy' . In the case of isolated peri- 



odic orbits the integrals over y and y' are evaluated by 
the stationary phase method while the integrals over z 
and z' are performed exactlyi&iii. This results in 



7 r— 1 



I'V V 2r e/ 

cos(rS 1 (Ep)/h— irrfi^/2) 



Cj(v x ,v x ) 



|det(M 7 -l)| 1 /2 



(10) 



where C 7 (i> x , v x ) is a correlation function of the longitu- 
dinal components of the velocity along the orbit 7 



Cj(v x ,v x ) 



dt 



dt' ^(iK(<+i')e^' /Tc! . (11) 



Formula i|l(J|l applies to the fully chaotic case where all 
orbits are unstable. In a lattice there are many copies 
of an orbit, so when counting different orbits one has to 
multiply each one with a degeneracy factor, depending 
on the geometry. Finite temperature and spin effects 
lead to additional factors of {irrT~/h0)J sinh(7rr T 7 /7i/3) 
and 2 cos(rT 7 ^si?/?i), respectively, where /is = eh/2m e 
is the Bohr magneton. There are similar formulas for the 
Hall conductivity g x; j 3 i 18 i 19 . 

Close to a bifurcation the saddle point approximation 
leading to (fTTTf) breaks down because the saddle point is 
not isolated. There are nearby saddle points from the 
other periodic orbits that participate in the bifurcation. 
This occurs when, by changing the energy or parame- 
ters of the system, the eigenvalues of M r — ► 1 and 11U|) 
diverges. In order to obtain the correct semiclassical con- 
tribution one has to integrate over the neighbouring sad- 
dle points as well. This results in transitional or uniform 
approximations for the bifurcating orbits. 

Bifurcations occur in specific forms that depend on the 
repetition number r for which det M r = 1 . They are 
characterized by normal forms which describe the charac- 
teristic motion of trajectories in the vicinity of a periodic 
orbit^. Before describing the normal forms we transform 
the integral in JSJ) into a form that is more appropriate 
for treating bifurcations. 

One of the integrals over y and y' is done in station- 
ary phase approximation, and the other is transformed 
to an integral over normal form coordinates^^. After- 
wards the integral over the z and z' coordinates can be 
performed. The resulting contribution from the r-th rep- 
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where 



(12) 



g ( E f) = jpf d ®' dp ex P (^( Q '' p )) • ( 13 ) 

Here So is the action of the periodic orbit at the centre 
of the Poincare surface, and v r is the number of conju- 
gate points along r repetitions of this orbit. 4>(Q'P) is a 
generating function for the Poincare map in normal form 
coordinates from (Q, P) to (Q' , P'). At the position of a 
periodic point 

-=P'-P = 0, ^ = Q-Q' = 0, (14) 

and $ is stationary. 

Eq. (|12f) is a transitional approximation for the bifur- 
cating orbits. It is correct when the orbits are close to 
a bifurcation, where its semiclassical effect is largest. 
If, due to a change of parameters, the orbits move fur- 
ther apart, it reduces to a sum of single terms as in 
(|10(l where, however, the amplitudes for the neighbour- 
ing orbits are inaccurate. The transitional approximation 
is sufficient for our purpose. A uniform approximation 
which has the correct single orbit limit, can be obtained 
by including the correlation function Cj(v Xl v x ) into the 
integral (where it is evaluated along neighbouring trajec- 
tories) and including also the Q' and P dependence of 
the pre-exponential factor that comes from the determi- 
nants In H12|) this factor has been evaluated at 
the central orbit and has the value one. 

Let us look at an example of a pitchfork bifurcation 
with normal form 



+ cC 



(15) 



where c > and the primes have been dropped for 
convenience. For e > there is only one solution of 
= = at (0, 0) corresponding to the central pe- 
riodic orbit. If, by changing parameters, e goes through 
zero, two new solutions appear which are located sym- 
metrically about Q = 0. In the generic situation (in 
systems without symmetries) this is a period doubling 
bifurcation. It occurs if the lowest repetition number r 
for which det(M r — 1) = is r = 2. At the bifurcation 
one new orbit of double the period arises which intersects 
the Poincare section twice. 

The integral (H2J can be evaluated in closed form for 
the normal form l|15l) and is given by a sum of two Bessel 
functions with index 4 and — j. At the bifurcation, e = 0, 
it has the simple form 



G[E) 



7rr(l/4)exp(37ri/4) 

ft 5 /4 c l/4 



(16) 
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The exponent of fi in the denominator, j3 = 5/4, shows 
that the bifurcation term is by a factor ftr 1 ^ stronger 
than that of a single orbit. The contribution of pitch- 
fork and tangent bifurcations to the conductivity were 
evaluated in£*2i. 

We note that the constant c in the normal form that 
appears in Ijlbj) can be obtained by following the actions 
and stabilities of the orbits through the bifurcations. It 
follows fromSi 



TrMo-TrMi 6e, S - S 1 



£ 

4c 



Hence one finds 



c = lim 



(TrMo - TrMi) 2 
144(^ - 5i) 



(17) 



(18) 



There is a second exponent, 7, which characterizes the 
semiclassical importance of bifurcations for spectral and 
transport properties. It specifies the region in parameter 
space (e in the example) in which the contribution of the 
bifurcation is strongest. Both [3 and 7 can be obtained 
from the normal form by a scaling argument. We perform 
a scaling of Q, P and e such that the integral in (|13|l takes 
the form 



G(E f ,e,h) = ±g(E f , 



1 



(19) 



In the example of the pitchfork bifurcation the argu- 
ment of the exponent is made fi-independent by changing 
the integration variables and a subsequent scaling of e. 



P 



Q 

-1/2 



1/4 



(20) 

5/4, as before, and 



As a result (3 = 
7 = 1/2. 

Different bifurcations have different exponents (3 and 
7 and, depending on the quantities that one considers, 
they are of different importance. We will consider in 
the following statistical properties of the conductivity. It 
will be shown that they are dominated by certain kinds of 
bifurcations, leading to universal properties in the regime 
of mixed classical dynamics. Before we do that we give 
an overview of classifications of periodic orbits that are 
known presently. 

The generic bifurcations of codimension K = 1 are 
those that occur if one parameter of the system is varied. 
They have been classified by Meyer and Bruno^S^. A 
list of the normal forms can be obtained from Table [I] 
by setting X2 = 1 (except for r = 1 where it is P 2 + 
X\Q + Q 3 ). These normal forms are simplified versions 
in which all constants and terms that are irrelevant for 
the determination of and 7 have been removed. The 
scaling proceedure of Eq. (|19fl yields the exponents j3 r 
and 7 r that are given in Table ITT1 

If more parameters than one are varied then other bi- 
furcations occur in which more complicated configura- 
tions of periodic orbits coalesce. For example, by vary- 
ing a second parameter one can make certain codimen- 
sion K = 1 bifurcations occur at the same instance. The 
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2 
3 
4 

5 
> 6 



P 2 + x x Q + x 2 Q 2 + Q 4 
P 2 + Xl Q 2 + x 2 Q 4 + Q 6 
(P 2 + Q 2 ) 2 + zi(P 2 + Q 2 ) + x 2 Re[(P + iQf] 

P 2 Q 2 + Xl {P 2 + Q 2 ) + x 2 (P 2 - Q 2 ) 2 
Re[(P + iQf] + xi (P 2 + Q 2 ) + x 2 (P 2 + Q 2 ) 2 
(P 2 + Q 2 ) 3 + xi(P 2 + Q 2 ) + x 2 {P 2 + Q 2 ) 2 



TABLE I: The relevant parts of the normal forms for bifur- 
cations of period-r orbits with codimension K = 2. 
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3 
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5 




1 
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2 


2 


5 
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2 


1 

2 


5 
3 


1 



TABLE II: Exponents Pk,t and 7if,r for generic bifurcations 
with codimension K = 1 or K = 2. 



codimension K corresponds to the number of parame- 
ters that are varied. Generic bifurcations of codimension 
K = 2 have been classified by Schomerus^ and the nor- 
mal forms are given in TableQ] For each parameter Xj one 
has an exponent Oj that specifies the range over which 
the bifurcation is important, Axj ~ h r7j . The total vol- 
ume in if-dimensional parameter space then scales as K 1 
where 



K 



(21) 



The exponents @K,r and jK,r f° r bifurcations of codi- 
mension K = 2 are also listed in Table ITT1 

Although bifurcations of higher codimension do not oc- 
cur, in general, when only one parameter is varied they 
still affect semiclassical approximations, because of their 
finite extention in parameter space. This is one of the 
main reasons why the semiclassical analysis of mixed sys- 
tems is so intricate. 

For bifurcations of codimension K > 3 there is no com- 
plete classification. Only partial results exist. Our main 
interest lies in bifurcations of arbitrary codimension K 
with repetition number r > 2K + 2. These have the 
(simplified) normal form 

K 

$ r , K (Q, P) = I K+1 + x « /n + 0( jK+2 )> ( 22 ) 

71=1 

where I = Q 2 + P 2 . Expressing $ in terms of Q and P, 
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we find 

0K,r = 2 



2 2K+1 

2(K+1) ~ K + l 
n 



K 



K + l 



sr K 

lK,r = 2_^i a K,r,n = — 
n=l 



(23) 



It will be shown that these bifurcations are the most im- 
portant for determining conductance fluctuations in the 
next section. 



III. MOMENTS OF THE CONDUCTIVITY 
FLUCTUATIONS 

We estimate the semiclassical size of the conductivity 
fluctuations a^ x by evaluating the ^-dependence of the 
moments 



M 2m (h) = ([a£(E,h)] 2m ) 



(24) 



where (•) denotes a local average over one or more pa- 
rameters of the Hamiltonian, for example Fermi energy, 
magnetic field or system parameters. In this section we 
assume that very long orbits do contribute to the longi- 
tudinal conductivity. The consequence of the damping of 
the contributions of long orbits due disorder and inelas- 
tic processes is discussed in the next section. The central 
point now is to replace the average by an average over 
the parameters in the normal forms 



M 2m , r , K (K) = B J d K x[C, 



fi 

r,K 



(x,fi)] s 



(25) 



where B is a normalization constant. The scaling pro- 
ceedure in (|19l) yields that each bifurcation, labelled by r 
and K contributes a term that scales as l/jT, 2m ^ K - r ^ K ' r . 
With the H dependence thus extracted, these contribu- 
tions can now be compared for different bifurcations. The 
bifurcation that wins the competition is that for which 
the h exponent in the denominator is largest, and it deter- 
mines that rate at which the x-averaged moments diverge 
in the semiclassical limit. That is 



M 2m (h) 



as 



0. 



where 



max(2m/3if r 

K,r 



7K,r) • 



(26) 



(27) 



The exponents v m are universal numbers that are de- 
termined by studying the hierarchy of bifurcations. Simi- 
lar universal exponents have also been found for moments 
of the fluctuating parts of the spectral counting function 
and wavefunctions2Si2SiSi. They have been named 'twin- 
kling exponents', in analogy to the exponents that control 
the intensity of twinkling starlight. 
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TABLE III: 


Values of v m 


,K,r — 


2mf3 K ,r 


- 7K,r- 
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1 2 


3 4 
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6 
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v m 




9 12| 


16 


19| 


23± 


26| 


K 


1 2 


2,3 3 


3,4 


4 


4 


5 



TABLE IV: Exponents v m and codimensions K of the domi- 
nanting bifurcations for generic two-dimensional systems. 



For the generic bifurcations with K = 1 and K = 2 
the exponents v m .K. r — 2m(3K,r — iK.r are listed in table 
TTT1 For generic bifurcations with K > 2 and r > IK + 2 
the exponents v m k r follow from Eq. (|23|l and are given 
by 



V m ,K,r = 2m- 



2K+ 1 
K+l 



K 



(28) 



which arc independent of r. Although the normal forms 
for K > 3 and 1 < r < K have not been classified com- 
pletely, it has been argued in refsi^SiSi that these bifurca- 
tions cannot contribute to the twinkling exponents since 
they have a counterpart with K < K and f > 2K + 2 
with the property that v m ^ - > v m ,K,r- The same ar- 
gument can be applied in the present case, and hence 
the maxima of the exponents given in l|28JI represent the 
maxima with respect to all generic bifurcations and are 
universal. They can be written as 



4m 



2m 

mm , 

K \K 



1 



K 
~2 



(29) 



The exponents v m are given in Table iTvl If m is a square 
then v m = 4m— 2^/rn + 1/2. If not, the maximum value 
is attained at one or both of the two integer values of K 
that are closest to 2y / m— 1. It is clear from the moments 
that the conductivity fluctuations do not have a Gaussian 
distribution. 

Another quantity explored in the literature concerns 
the magnetic fingerprint in antidot lattices incorporated 
in the auto correlation function of the conductivity 3 ^. 
Consider 



F(Ax) = (4 l x (x + Ax) x a£(x)) x 



(30) 



where a[ l x is written as a sum over classical orbits and x 
is a parameter of the system, e.g. the magnetic field. 
The average {')x.y might be over other parameters as 



6 



well. At Ax = this is identical to the second moment 
M 2 . The bifurcations then determine the characteristic 
length scale over which the correlations decay as Ax is 
increased. In the generic situation all parameters x n in 
the normal form are affected when the physical param- 
eter x is changed. Hence the correlation length scales 
as h a where a is the minimum of the o^. For generic 
bifurcations with r > 2K + 2 this yields a — 1/(K +1). 

IV. DISCUSSION 

The semiclassical analysis of the conductivity fluctu- 
ations leads to a powerlaw dependence of the moments 
M 2m ~ \/K Vm , with fractional exponents v m . In situa- 
tions where arbitrarily long periodic orbits contribute to 
the conductivity these exponents are universal numbers 
that are obtained from the competition between different 
bifurcation. The analysis was done in terms of H. If a 
different semiclassical variable is used instead of h, e.g. 
the Fermi wave length, then the exponents have to be ad- 
justed depending on how the classical action scales with 
the semiclassical parameter. 

In experiments most of the long orbits are suppressed 
and the conductivity is determined by a relatively small 
number of periodic orbits. One cannot expect then, in 
general, to see the universal exponents in particular for 
the higher moments, because they originate from high 
repetitions of periodic orbits. Instead the exponents 
will be determined by the competition within the much 
smaller class of those bifurcations that affect the relevant 
periodic orbits. The dominating exponent can then be 



found by comparing the corresponding exponents v m ,K,r 
in Table lllll (assuming that higher codimensions are not 
important). Nevertheless, also in this case the conductiv- 
ity fluctuations would be non Gaussian and the moments 
would have a powerlaw dependence with fractional expo- 
nents in clear contrast to the chaotic case with Gaussian 
fluctuations and v m = 2m. Also the length scale of the 
decay of correlations is larger if the classical dynamics is 
mixed, and is of order h a with a < 1. 

There is, however, one property of antidot lattices that 
might make it possible to observe the universal exponents 
for the lower moments, and that is their high symme- 
try. In systems with discrete symmetries there are also 
non-generic bifurcations which typically have the same 
normal form as generic bifurcations, but at lower repeti- 
tion numbers. For example, the pitchfork bifurcation Eq. 
(|I5|I can occur at repetition number r — 1 where two new 
orbits of the same length appear instead of one orbit of 
double the length 40 . Similarly, if the antidot lattice is 
e.g. invariant under rotations of 7r/3 then a bifurcation 
that is generic for r = 6 can occur at the first repetition 
of a periodic orbi^. This might make it feasible to ob- 
serve the universal regime that in non-symmetric systems 
is restricted to very long orbits. 
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